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Abstract 

The avalanche polynomial on a graph, introduced in [S], capture 
the distribution of avalanches in the abelian sandpile model. Studied 
on trees, this polynomial could be defined by simply considering the 
size of the subtrees of the original tree. In this article, we study some 
properties of this polynomial on plane trees. In [5], they show that two 
different trees could have the same avalanche polynomial. We show 
here that the problem of finding a tree with a prescribed polynomial is 
NP-complete. In a second part, we study the average and the variance 
of the avalanche distribution on trees and give a closed formula. 

1 Introduction 

Self-organized criticality is a concept introduced by Bak, Tang and Wiesen- 
feld [1] to describe the behavior of natural systems like earthquakes [H [16] , 
forest fires. A simple model that verifies this paradigm is the Abelian Sand- 
pile Model on the 2-dimensionaI lattice |11|, [TJ [10] . This model is based on 
a cellular automaton where each cell has a number of sand grains on it and 
one cell topples whenever the number of grains is greater or equal to four. 
In this toppling, the cell gives one grain to each of its neighbor. Thus, some 
other cells may topple and the sequence of toppling is called an avalanche. 

This model was also considered by combinatorists [21 [3l [6l [1^ but on 
general graphs. In |5j, a polynomial was introduced to encode the distribu- 
tion of avalanches that was previously studied from another point of view 
[T| [9l [8l [T3l I15j . In this paper, we study the inverse problem. Given the dis- 
tribution D of avalanches, can we find a tree whose avalanche distribution 
is D ? Moreover we study the average and the variance of this distribution 
on plane trees. 
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2 Avalanche polynomial on a plane tree 



2.1 Definition and main result 

The avalanche polynomial [5j encodes the size of avalanches in the abehan 
sandpile model. Let T = (V, E, r) be a tree rooted at r G y whose vertex- 
set IS V = {ui, . . . and edge-set E. Without loss of generality we 
can take r = Vn+i- Let v ^ V he a, vertex. A subtree T' rooted at u is a 
tree whose vertices are descendants of v in T. The maximal subtree is the 
subtree containing all descendants of v. The size of a tree is its number of 
vertices. 

We label the vertices Vi according to the following algorithm [5] (see 
Figured]): 

1. Label the root 

2. The label of a child v of a vertex labeled ^ is ^+|maximal subtree rooted at v\ 

In [5], the authors show that the avalanche polynomial is also defined by 
AvT{q) = 'Ylii>iP'i.^'' 1 where pi is the number of vertices labeled i in T. In 
the sequel we will either speak of the polynomial or its sequence of labels. 




AvT{q) = q"" + Iq"" + Zq' + 3g« + 2q^ + 4g^" + 

Figure 1: Example of an avalanche polynomial of a tree 

In [5], they show that the avalanche polynomial is not a tree invariant 
by exhibiting two different trees with the same avalanche polynomial. In 
the sequel, we study the problem of finding a tree for a given polyonomial. 

Theorem 1 Given a polynomial P{q), finding whether there exists a tree T 
whose avalanche polynomial is P{q) is NP-complete. 
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2.2 Proof of theorem [T] 

We show that the 3-PARTITION problem can be reduced to our problem. 
Definition 1 The 3-PARTITION problem In,c is the following: 

• INPUT ; C G N, (oi, 02, • • • , asn) be non-negative integers such that 

• OUTPUT: 

1. A partition of the integers into n parts of equal sum C. 

2. NO if there is no such partition. 

This is a well-known example of NP-complete problem ( [12J ) . Our proof 
has several steps: 

1. We associate a polynomial P to each instance of 3-PARTITION problem 
(see Definition [2]) . 

2. We show that P is an avalanche polynomial by exhibiting a tree Tp 
with AvTp (q) = P- Furthermore, the structure of Tp gives a solution 
to 3-PARTITION. (see Lemma [I]) 

3. Then, we prove that Tp is the unique tree whose avalanche polynomial 
is P (see Lemma El El d [5|) . 

Definition 2 Let In^c be an instance of the 3-PARTITION problem on oi, ... , 03^. 
P{q) is the polynomial associated to In,c' 

3n 3n, 

P{q) = nq^+^ + ^ q^+^+''^ + J^(ai - l)qC+-^+^ 

2=1 i=l 

Lemma 1 The polynomial P{q) associated to an instance In,c of the 3-PARTITION 
problem is an avalanche polynomial. Moreover, the tree given in Figure \M 
yields a solution to the 3-PARTITION problem. 

Proof 1 LetHk = {a{k), (3{k),j{k)) be a partition of {1, ... ,3n} in n parts 
such that 

O'a(k) + 0./3(k) + O'.y(k) = C 

Note that P{q) is the avalanche polynomial of the tree Tp in Figure [3 

The nodes at height 1 in Tp have 3 children, whose labels are a part of 
Hfc. Hence, there exists at least a solution to the Avalanche Polynomial 
Problem which yields a solution to the 3-PARTITION problem. 
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Figure 2: A solution to the 3-PARTITION problem 



Nevertheless, many solutions (different non-isomorphic trees) to the Ava- 
lanche Polynomial Problem might exist. 

To prove Theorem [H we ensure the unicity thanks to the following re- 
mark: each value Oj (and C) can be multiplied by a polynomial integer factor 
A = A(?7-). The avalanche polynomial becomes: 

3n 3n 

1=1 1=1 

Note that A is polynomial ensures that our transformation is a polynomial- 
time reduction. With the following substitutions (oj Xai and C — > AC), 
the tree of Figure [2] is a solution of P\. We are going to prove that, taking 
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A large enough, every tree such that Avriq) — ^xil) yields a solution to the 
3-PARTITION problem. 

Lemma 2 The root has n children, each one labeled AC + 1. 

Proof 2 By definition, the smallest value is the label of a child of the root. 
Hence then vertices labeled XC + 1 are children of the root. Butn{XC + l) = 
Aaj) +n, the number of vertices. Hence all other vertices are children of 
one of the vertices labeled AC + 1. 

Let / be the set of all vertices labeled AC + 1 + Aoj and L the set of 
those labeled AC + Aa, + 2, A > 2. Note that if A > 1 then \I\ = 3n, 
\L\ = n{XC -3) and I f]L = 0. 

Lemma 3 If X > 2 then the leaves of every solution tree are in L. 

Proof 3 Let v be a leaf of a solution tree labeled ly. Then its parent p is 
labeled Zt, — 1. If 1^ = AC + 1 + Xa^ there does not exist any vertex labeled 
AC + Xa^. Otherwise AC + Xa^ = AC + Aop + 1 + 6 (6 £ {0, 1}). Hence, 
X{ay — Op) = 1 + 6 which contradicts X > 2. p cannot be labeled AC + 1 since 

XOy 7^ 1. 

Lemma 4 // A > 3n then all vertices of L are leaves. 

Proof 4 Suppose that there exists an internal vertex in L. Take v as the 
one with the greatest label 1^ = AC + Xa^ + 2. Let u be a child of v. If v is a 
leaf, then its label is ly + 1 = AC + Xa^ + 3. But there is no such vertex since 
A > 3. Hence v is an internal node. But v is the greatest internal node of 
L and the label of v is greater than 1^, hence u £ I . 

Then v is either a child of the node labeled AC + 1 or not. The two cases 
are shown in Figure\^ 

Suppose that the parent p of v is labeled AC + 1 + Xop + (5, 5 = 0, 1. The 
size a of the subtree rooted at v is given: 

1. either by the difference between its label and the label of its parent, 

a = XC + Xay + 2- (AC + l + Xap + 5) = A(a^ - Op) + (1 - 6), 

2. or by the sum of the size of the t child-trees of v plus 1 (to take v into 
account): 

t t 
a = 1+Y^ ((AC + Xos, + 1) - (AC + Xay + 2)) = 1-t+Y^ A(o,,-a^) 

i=l i=l 
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Figure 3: Two different shapes of trees that appear in Lemma [H 



Hence A(a„ — a^) + {\ — 5) = 1 — t + Yl\=i ^i'^s^ — clv), and 
^ (ELi "s, - {t + l)av + ap) =t- 6. 

But t < |/| = 3n (the number of values ai). Hence t — 6 < 3n. Taking 
A > 3n concludes whenever p is labeled AC + Acp + 1 + 6,6 = 0,1. 

Note that when p is labeled AC + 1, taking Op = and 6 = in the 
preceding proof concludes. 

Lemma 5 // A > 3n, then the nodes of I are children of one of the vertices 
labeled AC + 1. 

Proof 5 Let v £ I be a vertex labeled AC + Xoy + 1. By Lemmas 13\ and^ 
the Xoy — 1 vertices labeled AC + Xa^ + 2 are leaves, hence children of v. 
Note that if k vertices vi, . . . ,Vk have the same label than v then there are 
k{Xay — 1) vertices labeled AC + Xa^ + 2. These vertices are children of the 
Vi,! < i < k. The size of each maximal subtree rooted at Vi is less or equal 
than Xa-u since the size of the subtree is the label of v minus the label of its 
parent p, and the label of p is at least AC + 1. Hence each Vi has exactly 
Xoy — 1 children labeled AC + Xoy + 2. 

Suppose that the parent p of v is labeled AC + 1 + Xap then the size of 
the subtree rooted at v is A(a„ — Op) < Aa„. 

Thus the parent of v is labeled AC + 1. 
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2.3 Avalanche polynomial on plane tree of height at most 2. 

We showed in the last section that finding a tree whose avalanche polynomial 
is given is NP-complete. In this part, we show that considering only trees of 
height at most 2, then the problem becomes poynomial. Moreover we give 
a linear time and space algorithm to solve it. 

A tree of height 2 and its labeling is given in Figure U] and its associated 
polynomial is P{q) = XlILi i^""' + i^i ~ ^)q"''^^) ■ Note that if one of the 
is equal to 1 then it is a leaf rooted at the vertex labeled 0. 

The following algorithm takes as input a polynomial P{q) = Y17=i 
and output either: 

• A tree T of height at most 2 whose avalanche polynomial is P. 

• NO if no such tree exists. 

1. Create the root with label 0. 

2. Find a child of the root by taking the first non-nul coefficient aj. If 
j = 1 we have found aj leaves of the root. If j > 1 we have aj nodes. 
Each of these nodes must have j — 1 children labeled j + 1. Hence flj+i 
must be greater than aj{j — 1). If not, output NO and exit. 

3. Subtract ajq^ + aj{j — l)q-^^^ from the polynomial; we take this new 
polynomial as P 

4. If P = output the tree and exit. 

5. Goto El 




Figure 4: Labeled tree of height 2 
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3 Avalanche polynomial on plane trees 

In this section, we study the avalanche polynomial on plane trees. The aim of 
this section is to find the average distribution of avalanches on a random tree. 
We first give the closed expression for the number of avalanches of a given 
size in plane trees. Unfortunately, this formula does not help to retrieve 
informations on the distribution, thus we give a recursive formula for the 
distribution and deduce the mean and the variance of this distribution. 

3.1 Definition 

Let us note Ck = i^j^) the A;-th Catalan number. Let 7^ be the set of 
rooted plane trees with n edges. |7^| = C7„. We extend the definition of 
avalanche polynomial to the set of plane trees (see Figure [5]) : 

^(*' <?) = E E ^^^(^)*^ = E M^)^' (1) 

j>0 TeTi p>0 



Size 


Trees 


Polynomial 


1 








2 




qt 
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{2q + q^ + q^)t^ 
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{5q + 2g^ + 4q^ 
+2q'^ + q^ + 



Figure 5: First terms of A(t, q) 



Remark 1 Note that : 
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In fact, the coefficient of in A{t, 1) is the number of rooted plane trees 
of size p+1 times the number of vertices (excepting the roots). Hence it is 
pCp. 

3.2 Closed formula 

In this section, we exhibit a decomposition of the rooted plane trees which 
yields a closed formula for the avalanche polynomial. 

Theorem 2 The coefficients of the polynomials An {q) are given by: 

L 2 J / k \ 

mAn{q)= J2 E Cp,.AiiCp^.p,_Acn-p,+i (2) 

k=i pen„^„^fc \i=2 / 

where 

k 

^,v,k = {{Pi,---,Pky,Pi <Pi+i,Pk < n,^pi = v} 

i=l 

Proof 6 The coefficient [q^]An{q) represents the number of vertices labeled 
V in every plane trees of size n + 1 . 

Let X be such a vertex labeled v in a tree T and yo = x,yi, . . . ,yk = Xn+i 
be the path joining the root Xn+i labeled to the vertex x in T. Let lab be the 
function which maps any vertex y ofT onto lab{y), the label of y iuT. The 
proof is based on the study of rooted plane trees which have a given path. 



• Note that k < 



-l+Vl+Sv 
2 



because the labels of the vertices of the path 
from the root to x take minimal values when the maximal subtree of T 
rooted at yk-i is reduced to the path yo, . . . , yt-i- 

Then the labels of the vertices on this path are k, k + (k — 1) ,k + (k — 
l) + {k- 2), ^^^^ and if k > -^+^+^ then the last label is 
stricly greater than v. 



The sequence of labels of the vertices yi is stricly decreasing by defini- 
tion. 

The sequence of the differences {lab{yi) — lab{yi^i)) is stricly increas- 
ing. This relies on the fact that each difference is the size of the maxi- 
mal subtree rooted at yi^i and that the subtrees are nested. Moreover, 
lab{yk-i) — lab{yk) < n as this is the size of the maximal subtree rooted 
at yk-i- Finally, Yli=i ^a^(yi-i) - laKVi) = lab^yo) - lab{yk) = v. 
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Let k be an integer such that k < i+\A+S '£ [pi, . . . ,pk) be 

a stricly increasing sequence such that Pk ^ n and Yl\=iPi ~ ^- ^ 
be a path of vertices {yo, - ■ ■ ,yk} labeled lab{yo) = v,lab{yi) = lab{yo) — 
Pi,lab{y2) = lab{yi) - p2, . . . ,lab{yk) = lab{yk-i) - Pk = (see Figure\B^. 
Denote by Fi^i (resp. Fi^2) ths forest hanging on yi at the left (resp. right) 
of the edge {yi,yi-i). 

Lemma 6 The number of plane trees of size n + 1 having P as subpath is 

Cpi-l (Y[i=2 ^Pi-Pi-i) C'n-pfc+l- 

Proof 7 Note first that the size of the maximal subtree hanging on yi is 
Ui — Vi+i = Pi+1- Then the number of vertices of the forest Fi^i U Fi^2 is 
Pi+i — Pi — 1 for i € {1, . . . , — 1}. The number of such pair of forests 

{-^i,i,-^i,2} is C'pi+i-pi- 

The size of the maximal subtree rooted at yg is yo ~ Ui = Pi- There are 
Cp-^^i such trees. 

The remaining nodes are in the forests F^^i and Fk^2 which union has 
n — Pk vertices. So that the number of such forests is Cn-p,.+i. 




Figure 6: A tree with its subpath P 



Remark 2 The coefficient oft^q in A{t,q) is Cn and the coefficient oft'^q^ 
is Cn-l. 

Proof 8 There are two different proofs of these results, the first one consists 
in taking v = \ or v = 2 in Equation 

The second one relies on the following property: 

Lemma 7 Let T' be a tree. Let S^' be the set of trees T of size n + 1 where 
T' is a subtree of T rooted at the root of T . Then \St>\ = Cn-\T'\+2- 
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Figure 7: Bijection between St' and plane trees 



Proof 9 The bijection between the trees of St' and plane trees of size n — 
|r'| + 1 is given in Figure^ 

For V = 1 take T' as a single edge. For v = 2 take T' as a path of length 

2. 

3.3 Functional equation 

Lemma 8 The polynomial Ap(q) , p > - the distribution of avalanches on 
plane trees of size p+ 1- is defined by the following recurrence: 

p 
fc=0 

This yields a functionnal equation for A{t, q): 

Proof 10 The proof relies on the following decomposition of rooted plane 
trees (see Figure IB^. A tree Tp+i of size p + 2 can be decomposed into a 
tree of size k + 1 and a tree Tp^k of size p + 1 — k as shown in Figure 
[3 Hence, the root of the tree Tk is the leftmost child of the root ofTp-.^ in 
Tp+i . Furthermore, the root of Tk is labeled k + 1 the number of vertices of 
Tk. This vertex contributes to the term CkCp-kQ^^^ in Ap^i{q). 

The avalanche polynomial of Tk is Ak{q). In Tp+i the label of a ver- 
tex of the subtree Tk is increased by k + 1. They contribute to the term 
Cp-kQ^^'^^^kiQ) ^p+i{q)- I'he last term takes into account the labels of 
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the vertices ofTp^k which stay unchanged since the root of Tp^k is the root 
of Tp+i . 




Figure 8: Decomposition of a plane tree 



Substituting Equation ^ in Equation ([T]j yields a functionnal equation 
forA{t,q): 



A{t,q) 



(1 - Vr^)(l - VI - 4gt) 



At 



p>0 fc=0 

(1 - \/r^)(i - VI - 4gt) 



+1 



4t 

1 - VT^^ 



, 1 - VI - 4t 



2 ^' 2 

3.4 Average of the distribution of avalanches 

Proposition 1 The average size of an avalanche on a plane tree of size n 
is asymptotically equal to ^r?!'^ . 

Proof 11 

n ra 

The mean size of avalanches on all plane trees of size n is given by 



M„ 



[Mt,l)]n 



TlCfi 



(5) 
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Equation ^ yields: 

Factorization o/^(t, 1) anrf substitution of A(t,l) and ^(t, 1) wzi/i 
t/ieir values of Remark [7] yield: 



dA ^ (7r^-i)(2t-i) 

5g ^ ' ^ 2(4t - 1)2 

Since 

j>0 



|^(t,l) = tj;4Hi + l)t'+ J](C7,+i-2Q)4^(j + lK+^+2 

i>0 i,j>0 

= tY,^\^ + l)f + Y,^H'^'E(a^^-^c^)w^ 

i>0 k>l i=0 

^ Cj+i - 2Ci _ kCk ^ i(C,+i - 2Ci) _ (2 + 5fc + fc2)Cfc ^ 

3 A - 

— (t,l) = ^(4*~i(f + 2)-(-l + i + 2i2)a-i)i^ 

^ i>i 

Substitution o/^(t, 1) in Equation ^ yields: 

^ 4»-Hn + 2) _ (n + l)2 
nC,i 2n 

Since Cn ~ M„ ~ ^n^/"^. 
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3.5 Variance 



The variance of the distribution is given by 



aq- 



nCn 



+ Mn- M?; where 



M„ denotes the mean size (computed in the last section). 



Hence we only need the expression of 



(1 - 4t)3/2 



(t,l) + 2-^(t,l) + 2t— — 



dq 



dtdq 



Thus we need the other derivatives ^{t, q),^{t, q), ^^{t, 1) and ^^(i, !)• 



Substituting their expression in ^^{t, 1) yields: 



dtdq 



d^A 



dq^ 
Since : 



(i,l) = -l/2 



(-1 + VI -4t) (10 i + Vl -4t - 1 - 32 12 + 32 13 
(-l + 4t)^ 



4" 



(1 - 4t)-4 = ^ (n + l)(n + 2)(n + 3)t" 

n>0 

n _ 4.^-7/2 _ V- (n + l)(2n + l)(2n + 3)(2n + 5)4"C„ ^, 
We obtain: 



n>0 



Vn 



73 



26 



— n H n H n -\ ^ 

15 60 15 15 2n 2n 
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n-l 



+ 



n2(n + l)2C2 
4"-i 



(n^ + 6n3 + ISn^ + 12n + 4) 



n2(n + l)C„ 
The leading term is: 



{n'^ + + 5n + 2) 
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n-l 



-n 



4 M 

n ~ - 



vr 



15 n2(n + 1)2(72 16 



n 



In fact the distribution show strange peaks and after a renormalization 
(scaling in x by a factor 1/n for the distribution on trees of size n and rescal- 
ing the ordinate between and 1) it gives the curves in Figure and EJ In the 
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following figures, we draw a curve for AiQ{q), A2o{q),.-. Moreover for such 
a polynomial An{q) = J2i>oPi'l^ '^^ P^^ points at coordinates {i/n,pi/pi) 
since pi is the maximum. Let's recall from Remark [2] that pi = Cn- The 
grey scale goes from white (for the first polynomial) to black (for the last 
one). 

I.On 1.0-1 




Figure 9: Distribution for the first 500 and 1000 polynomials 



4 Conclusion 

The Figures ([9]) and ([9]) point out some convergence of the renormalized 
distribution. A question arises naturally: 

• Is it possible to retrieve the asymptotic of the coefficients from their 
closed form ? 

Note first that the peaks that appear at integer abscissa x, x <C \/n do not 
converge to 0. In fact consider the trees made of a path ttq = 0, vri, . . . , vr^. 
of length X and a plane tree of size n — x rooted in tTx- The label of tt^ is 
n+(n — l) + (n — 2) + ...+(n — x+1) thus tends to x after the renormalization 
by n. Moreover there are Cn-x-i such trees. Since pi = Cn the peak is 
higher than Cn-x-i/Cn ~ (l/4)-i. 
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